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Abstract 

We study inclusive Vf, and cross sections in ^^C in a theory that takes 
into account significant nuclear renormalization of strengths. Our calculation is 
in excellent agreement with the measured inclusive muon capture rate and the 
flux-averaged Ve cross section, but overestimates the flux-averaged inclusive 
cross section. These reactions are of crucial relevance to the issue of possible 
to Ve oscillations. We also calculate the flux-averaged cross sections in ^^C and 
^^^Z, which are found to be consistent with the available experimental result. 
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I.INTRODUCTION 



Extracting an interesting physics signal of some rare events from an experiment 
crucially depends on our ability to understand the physics background in that setting. 
A good example of this is the search for neutrino oscillations, and hence physics be- 
yond the standard model, in the Liquid Scintillator Neutrino Detector (LSND) related 
experiments, with profound implications for particle physics, nuclear physics and as- 
trophysics. The LSND group have advanced evidence for the —>■ Ve oscillations 
in a recent experiment with neutrinos from pion decays in flight 1^. Our trusts in 
this claim depend crucially upon the performance of these experiments in benchmark 
reactions on nucleons and nuclei initiated by neutrinos. Some of these benchmark 
reactions are exclusive ones on proton, ^^C and targets, while others are inclusive 
processes in complex nuclei. Given separate claims of agreements with experiment 
and lack of it in recent theoretical calculations ^, ^, |^, ^ of inclusive neutrino 
reactions in ^^C, this is of topical interest and deserves a careful examination. 

In this paper, we are going to discuss inclusive neutrino reactions in nuclear targets 
^^C, ^^C and ^''A/, concentrating on the theory of the reactions |^ 

^'C + z/e ^X + e-, (1) 
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C + z/^ ^ X + ^" 



(2) 



where the produced X is not observed. Thus, generically, we can indicate the reactions 
(1,2) as ^'^C{i'i,l~)X where X is the unspecified nuclear state and / is the charged 
lepton. At low energies, the reaction (1) has been measured by the E225 and LSND 
1^, |TU| experiments at the Los Alamos Meson Physics Facility (LAMPF) and the 



experiments by KARMEN [|Tl|, Q collaboration at the ISIS facility, with t'gbeam from 
muon decays at rest (DAR), given by the Michel spectrum. The flux-averaged cross 
sections from the various experiments are given in Table I. Particularly interesting to 
us are the most recent values of the cross section for the reaction (1) 



a{ue) = (14.8 ± 0.7 ± 1.1) X 10"^^ cm^ 



(3a) 



from the LSND experiment fiy] and 



ally. 



(15.2 ± 1.0 ±1.3) X 10"^^ cm^ {3b) 

from the KARMEN collaboration |jl2|. In Eq.(3), the inclusive cross sections are 
obtained by adding the cross sections leading to ground state and excited states, 
with their errors added in quadrature. Within their errors, these experimental results 
overlap. The reaction (2), on the other hand, has been measured over the last five 
years at LAMPF |IB| , using ly^ beam from pion decays in flight (DIE). The 



earlier E764 experiment of Koetke et al. [|T3[ used a neutrino beam of slightly higher 
energy than used in later experiments and gave a cross section for the reaction (2) 
too large when compared to the theoretical predictions [0, H, ^, ^, 0. This reaction 
was further studied by Albert et al. ||T^ with more massive detectors and larger 
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exposure than used in [13|, using a beam of neutrinos with energy around 180 MeV. 



These authors obtained considerably smaller cross section than that by Koetke et al. 



The recent studies by the LSND collaboration |]15[ have improved the experimental 
situation, providing the latest measured value for this cross section. We give, in 
Table I, the results of all of these experiments, but use in our theoretical discussion 



the latest result, reported by the LSND collaboration [ff^ 



a{v^) = (11.2 ± 0.3 ± 1.8) X 10"^° cm^. (4) 

We also calculate the flux-averaged neutrino cross-sections in ^^C and '^'^Al as bench- 
marks for our theory. These are found to be reasonably consistent with the available 
experimental results |^. 

Current generation of calculations of the above reactions can be grouped into two 
classes, depending on their predictions for the fiux-averaged cross sections for the 
process (2): those that substantially exceed this observed cross section 0, ^, |^ and 
others that find agreement [|, ||. The deficit of the fiux-averaged cross sections may 
be: (1) a manifestation of theoretical problem, of not being able to do a correct enough 
nuclear calculation; (2) an experimental problem of not doing a precise and reliable 
enough experiment. Since both theoretical calculations and experimental analysis 
are involved in the determination of the excess events in the channel in the LSND 
neutrino oscillation experiment it is important to have a clear understanding of 
the nuclear physics related uncertainties in this reaction. The purpose of this paper 
is to narrow down the options by examining the first point very critically from our 
point of view . 

A reaction which can be regarded as a benchmark in the context of the processes 
(1) and (2) in general, and (2) in particular, is the nuclear capture of muons (NMC) 



from the atomic 15" state by the charged weak current |T^, p!7| : 

^^C + iJ-ilS)^X + u^. (5) 
This process serves as an excellent check, in the low and intermediate energy transfer 



region, in our ability to control the theoretical uncertainties [17|. The inclusive NMC 



rate Ac for the process (5) is very accurately known. Taking the world average of the 



best experimental determinations of the inclusive muon capture rate Ae,||T8|, |T9[with 
their errors added in quadrature, we obtain[0, 

Ac = (3.80 ±0.10) X 10^"^ (6) 

Thus, we have here a weak reaction rate, which is closely related to the processes (1) 
and (2) and is known at an accuracy of about 2.5 %, posing a tremendous challenge 
to the theoretical approaches to weak nuclear reactions in nuclei. 

Various theoretical approaches to calculate reactions (1) and (2) 0, ^, ||, |^, when 
applied to the inclusive muon capture, reproduce the NMC rate Ac quoted in Eq.(6) 



rather well, within the limits of their theoretical accuracy [T^, ^ . In the case of 
inclusive neutrino reactions, however, the situation is different. The calculations of 
Kolbe et al. reproduce the (fg, e~) rather well, but overestimate the (z/^, fi~) by about 
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50%. Auerbach et ai, on the other hand, use the set of parameters for their model, 
which explain the inclusive muon capture rate. They predict cross sections for the 
neutrino reactions, which are 20% and 15% larger than the experimental values for the 
reactions (1) and (2) respectively, quoted in Eqs.(3) and (4). These overestimates are 
from the maximum experimental values allowed within the quoted errors. Therefore, 
they represent a real discrepancy with the experiment. It is possible to explain the 
observed neutrino cross sections in the calculations of Auerbach et al. with another 
version of their residual nuclear forces by varying the model parameters, but this 
version predicts a NMC capture rate of 3.09 x 10^s~^, which is rather small compared 
to the value quoted in Eq.(6). 

Finally, to complete this survey, we would like to mention that a Fermi gas model 
calculation with a Fermi momentum kp = 225 MeV gives a much higher value of 
24.1 X 10~^° cm^ for the flux-averaged (f^, cross section. This is reduced to 22.7 x 
10"'^'' cm?, when the effects of meson exchange currents are taken into account 0. In 
another calculation , the method of so-called "elementary particle model", extended 
to inclusive reactions, has been used, to obtain a cross section of 13.1 x 10"^'' err? 0, 
which is in good agreement with the experimental value. However, the extension of the 
elementary particle model to the inclusive reactions makes use of several assumptions, 
which have not been tested in the energy region of the LSND experiments. Here this 



method is expected to underestimate the inclusive cross sections g |1^, p2[ . 

This survey brings us to our calculation, which we briefly describe in section II 
and discuss the results and conclusions in section HI. 

II. FORMALISM 

The matrix element for the neutrino nucleon reaction for a neutrino of flavor 
/, (/ = e, /i), i e. 



is given by 



where 



i^i{k) + n{p)^l-{k')+p{p'), (7) 



T=—j^u{k'mi-l,)u{k)J,, (8) 



J, = n{p')[Fn<f)l, + Fn^')^^,^f^ + FW)w + F^{q')%lMp)- (9) 

In Eq. (9), q^ = k^ — k'^, is the four momentum transfer, , F2 , F\ and F^ are the 
known weak nucleon form factors. The double differential cross section cro(g^,A;'), is 
then given by 

k' 
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where SS|Tp represents the sum and the average respectively over the final and the 
initial spins of the leptons and the nucleons and is evaluated exactly, using the matrix 
element T defined in Eq.(8). Its analytic expression is given in ref. 0]. In a nucleus, 
the neutrino scatters from a neutron moving in the finite nucleus of neutron density 
Pn{r), with a local neutron occupation number nn{p, t). Then the cross section in the 
local density approximation and in the free nucleon picture is given by 

cr(g^ k') = 2 J drj^^Hnip, r)cro(g^ k'), (11) 

where the neutron energy En and the proton energy Ep in the expression of (Jo{(fi k'), 
given in Eq.(lO), are replaced by Enijp) and Ep{p') respectively, p and p' being the 
momenta of the neutron and proton in the nucleus. However, neutrons and protons are 
not free and their momenta are constrained to satisfy the Pauli principle i.e.p < ppn 
and p' > pFpi where ppn and ppp are the local Fermi momenta, given by 

ppn{r) = [3^Vn(r)]^/' andp^p(r) = [?>7i^ pn{r)Y'\ (12) 

Moreover, in the finite nucleus, there is a threshold energy for the reaction to proceed, 
also called the Q- Value, and this should be taken into account. Finally, the charged 
lepton produced in reactions (1) and (2) moves in the nucleus and its energy is 
modified by the Coulomb interaction, which should be accounted for. In our approach, 
these effects are incorporated by modifiying the argument of the 5 function in Eq. 
(10), from E - E' + En- Epio E - E' - V;(f) + En{p) - Ep{p'), and replacing the 
factoi J j^nnip,r)^f^S{E - E' + En - Ep) occurring in Eq. (11) by -^ImU{qo,q), 
where 

qo = E-E'-V,-Q + Q'. (13) 

In Eq.(13), Vc{r) is the Coulomb energy of the lepton and Q' = Epp — Epn, is 
introduced to take into account the unequal Fermi sea in the case of N Z nuclei. 
U{qo, q, is the Lindhard function given by 

TT(n rr)= ( "^'^ nn{p)[l-np{p + q)\ 

J {2^f qo + En{p)-Ep{q + p) + ieEnEp ^ ^ 

With these modifications, the total cross section a{Ey) is given as 0: 
aiE,) = / r^dr / k'^dk' / (i(cos^)— — 

TT JO JpJ"™ J-1 EijEl 

xj:T.\T\HmU[E, -Ei-Q + Q'- K(r), ^Q[Ei + V,{r) - rm]. (15) 

The kinematic limits for the lepton momentum k' are easily computed in our 

special case|^]. For the numerical integrations, we use Gaussian quadrature with high 
enough accuracy for our purpose. The radial integration in Eq. (15) is performed up 
to a radius R = ci + 5fm, where ci, is the radius parameter in the two-parameter 
harmonic oscillator and Fermi density distributions used for the nuclei, considered in 
section III. 
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The renormalization of weak currents in the nuclear medium is taken into account 
by calculating the effect of propagation of the particle-hole (ph) excitations in the 
nuclear medium on various terms occurring in SS|Tp. The ph response is then re- 
placed by a Random Phase Approximation (RPA) response accounting for the ph and 
the A/i components, which interact through an effective spin-isospin nuclear interac- 
tion described by the Landau-Migdal potential. The details of this renormalization 
procedure as well as those of Eq. (15) are given in |^. However, we have made here 
the following improvements, which considerably reduce the theoretical uncertainties 
in our calculations from the previous versions of our model: (1) Our new Lindhard 
function makes use of a strategy |^ that avoids the pathologies of the ordinary 



Lindhard functions in the limit of go? ^ both going to zero, q'^ being the momentum 
transferred to the nucleus in the processes of interest. (2) The nuclear response func- 
tion is renormalized by the ph and the correlations in nuclei |]^, effects of which 
are quite large for low and intermediate energy neutrino scattering. The physics of 
this renormalization depends, among other things, on the Landau-Migdal spin-isospin 



parameter g' pj]. This itself has an uncertainty of ±0.1 around its preferred value of 
0.7 P5|. We take into account the theoretical uncertainties of our estimates of (j{Vf.) 
and cr(zy^) due to this variation of g' . (3) Finally, the target nucleus ^^C has intrinsic 
parametric uncertainties in the radial density function. The effect due to this uncer- 
tainty in our cross-section estimate is taken into account by repeating our calculation 



in several radial parametric settings p6[. Overall, we achieve a theoretical accuracy 
around ±10% for a{V(.) and cr(zy^). This significant improvement in theoretical ac- 
curacy is even better in the case of the NMC rate for the inclusive process (5) [|18 



We estimate here an uncertainty of ±2% due to nuclear radial effects and ±5% due 
to the variation of the spin-isospin parameter g' . Treating these two uncertainties 
independently and adding them in quadrature, we get a theoretical error of about 
±6% and obtain 



Ac = (3.60 ± 0.22) X 10V\ (16) 

it in excellent agreement with the precise experimental data (Eq.6). The inclusive 
nuclear muon capture provides us with a critical benchmark, an independent accurate 
check of our ability to describe nuclear inclusive weak processes clearly related to the 
neutrino scattering. 

In summary, our method used in this paper is essentially an RPA approach built 
up from single particle states of an uncorrelated local Fermi sea. This method is, in 
practice, found to be a very accurate tool, when the excitation energy is sufficiently 
large such that relatively many states contribute to the process, in particular, if a 
large fraction of it comes from excitation to the continuum, as it is in the present case. 
The adaptation of this method to finite nuclei via the local density approximation 
has proved to be rather advantegeous to deal with inclusive reactions and has been 
successfully applied to the photonuclear reactions [0, electron scatter ing|^8[, deep 
inelastic scattering pOf and muon capture |jl8|, ^ . 



The numerical evaluation of the neutrino-nucleus reaction cross section is done 
using Eq.(15) and the rsults are presented in section III below. 
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III. RESULTS AND CONCLUSIONS 



In order to compare with the experimental results of KARMEN||12|| and LSND |p 
collaborations, we compute the flux-averaged cross section, 

_ _ S^^^^ a{E,) uj{E,) dE, 
/I; u:{E,)dE, 

where the neutrino profile function uj{Ei^) is well-known (i.e., the Michel spectrum in 
Ei,^ and the spectrum of E,^^ provided by the LSND experiment). The lower limit Eq 
in Eq.(17) is taken to be zero for the [ug, e~) reaction and 123.1 MeV for the {ly^, /i^) 
reaction |ll3| . 

For ^^C, we present our results in Tables II and III, and compare with experiments 
and other recent theoretical works in Table IV. Here are the main points of our 
analysis. In both Tables II and III, the rows 1 through 4 indicate four different 
choices of the radial parameters for the nuclear density p6[. The columns represent 
different choices of the Landau-Migdal spin-isospin parameter g'. In Table III, we have 
used the i^^ spectrum reported in the LSND papers |[l|, |1^ for a direct comparison 
with the experiments. The radial uncertainties are typically about 2%, while the g' 
variation represents a ±7% spread, for the {ue,e~) case and ±8% for the {v^,^^) 
case both around the central values corresponding to g' = 0.7. Thus, the overall 
spread from the theoretical error, taking both of these effects in quadrature, is ±7.3% 
for the (z^ejC") case and ±8.2% for the {ly^, case. In Table IV, we compare the 
presently available theoretical results with the most recent experimental results for 
these reactions. 

From Table IV, we can provide our best estimate of a{ue) and o"(f^) as follows: 

a(z/e) = (15.48 ± 1.13) x 10"^^ cm^ (18) 



i-i-J 



(16.65 ± 1.37) X 10"^°cm^ (19) 

We are in excellent agreement with the experimental value of cr(z/e), but our lower 
limit cr(z/^) is 15% higher than the higher limit o/cr(z/^) measured by the LSND col- 
laboration. 

For ^^C, we obtain a flux averaged cross section of 7.25 x 10~^^cm^ for Michel 



spectrum, using the density distribution parameters given in 
pared with the calculations of Arafune et al. 
9.58 X 10 



31 



.This should be com- 
who obtain a cross section of 



-41 



cm 



for the transition to the ground state and first excited state of the 



final nucleus, which together give 85 ' 
elusive cross section of 11.3 x 10~^^' 



% of the total coss section. This implies an in- 
"cm^. There could be a reduction of (10 — 15) %,if 
the momentum dependence of the form factors are taken into account |22|. This value 
seems to be in agreement with an unpublished result of Donnelly quoted by Krakauer 
it et al. p[. The calcultions of Arafune it et al. do not take into account the 
possible quenching of the weak interaction operators in nuclei, which is studied by 
Fukugita it et al. |BB 



In this paper, the quenching of the matrix elements of the 
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Gamow- Teller (G-T) operators is obtained in an effective operator approach, which 
takes into account the effect of core polarization, isobar and the meson exchange cur- 
rent processes. This leads to a 20 % reduction in the flux-averaged cross section for 
the ground state transition, while the cross section to the first excited state is reduced 
by a factor 3. Assuming, as before, that these two states together contribute 85 % 
of the total cross section, a flux-averaged cross section of 5.4 x 10~^^ cm^ is inferred 



from the calculations of Fukugita et al. [Q. 

We find that our flux- averaged cross section for ^^C, reported in this paper, is 
35 % smaller than that obtained by Arafune et al. while it is 35 % larger than the 
results of Fukugita et al. It will be interesting to test these predictions by measuring 
this cross section in the low energy neutrino experiments with liquid scintillators, 
where ^^C forms part of natural carbon. In the experiments of Krakauer et al. P], it 
is reported that 

Gav = 0.723a(z/]^C) + a(z/f < 18.3 x IQ-^^cm^. (20) 
In our approach, we calculate the flux averaged cross section in '^'^Al to be 11.48 x 



10~ cm with g' = 0.7 and using one set of parameters from ref. ||26[. This, along 
with the value obtained for the neutrino cross section in ^^C, gives a value of 16.5 x 
10~*^ cm^ for aav in Eq.(20). We associate a theoretical uncertainty of 6 % due to g' 
and density variation of ^''A/ on this average cross section. This value is consistent 
with the available experimental information on these reactions. 

We would like to emphasize that the renormalization of nuclear strengths in our 
model produces a reduction of about 40% in the (i/e,e~) cross section to bring it 
in agreement with the experiment [Q. Similarly large reductions also occur in the 
(z/^, case. 

In summary, our calculations show no discrepancy with the measured flux-averaged 
Uf, cross sections in ^^C, like other authors. We also nicely reproduce the measured 
inclusive muon capture rate, now known very accurately [1^, But we see a discrep- 
ancy, at least by 15 % in the flux-averaged cross section compared with the LSND 



experiment |T3[, the theoretical prediction being higher than the experiment. The 
discrepancies between the experimental and various theoretical results for (z/^,/i~) 
inclusive cross sections in ^^C should be taken seriously, in view of its implications 
in present studies of neutrino oscillations. Our results for the case of ^^C and '^'^Al 
are consistent with the only experimental limit available at present. A high-quality 
experimental measurement of the inclusive cross section, specially in ^^C, will be very 
useful in understanding the quenching of the G-T strengths in this nucleus in the light 
of the wide range of theoretical predictions for this reaction. 
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Table I: A summary of flux-averaged and z/^ cross sections as obtained in various 
experiments done with ^^C target. The unit for the Vf. cross section is 10~^^cm^ and 



for cross section is 10 



0-(z/e) 


LSND collab. 
14.8±.7±1.1 [g 0] 


LAMPF E225 

i4.i±2.3 g 


KARMEN coll 
16.8±1.4±1.7 1 
15.2±1.0±1.3 1 


ab. 

12 




8.3±.7±1.6 |T^ 
11.2±.3±1.8 |15| 


LAMPF E764 
159±26±37 ig 





Table II: Flux-averaged z/g. Four radial parameter sets are chosen from the lit- 
erature with parameters Ci and C2 in fms. The Landau- Migdal parameter g' is 
taken as g' = 0.7 ± 0.1. The unit for cross section is 10~^^ cm^. 



Cl 


C2 


g' = 0.6 


g' = 0.7 


g' = 0.8 


1.687 


1.067 


16.87 


15.69 


14.72 


1.672 


1.150 


16.65 


15.49 


14.53 


1.649 


1.247 


16.30 


15.15 


14.20 


1.692 


1.082 


16.99 


15.80 


14.52 



Table III. Flux-averaged cross section. Radial parameter sets and g' values are as in 



Table II. The flux is taken from Si 



Tsi The unit for cross section is 10 '^^ 



cm 



Cl 


C2 


g' = 0.6 


g' = 0.7 


g' = 0.8 


1.687 


1.067 


18.35 


16.82 


15.61 


1.672 


1.150 


18.19 


16.67 


15.45 


1.649 


1.247 


17.91 


16.38 


15.17 


1.692 


1.082 


18.45 


16.92 


15.70 



Table IV. Summary of flux-averaged cross sections. Experimental results are inferred 
by adding ground state and excited state contributions for {ue,e~). Theoretical results are 
from Kolbe et al. Auerbach et al. Q, Umino et al. [^] and this work. The units are 



10- 



"^^cm^ for 



and 10 



^^CTTT? for 



fu,/U ) cross sections. 





Kolbe 
et al. g 

15.6 


Auerbach 
et al. § 

12.9 - 22.7 


Umino 
et al. [| 


This work 
15.48 ± 1.13 


Exp 

14.8 ± 1.0 ± 1.5 1 
15.2 ± 1.4 ± 1.8 


, 10| 
12| 




19.3 - 20.3 


13.5 - 15.2 


22.7 - 24.1 


16.65 ± 1.37 


11.20 ± .3 ± 1.8| 


15|| 
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